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ABSTRACT • . 

The analysis In this paper Is concerned with the plezovl scous-rlgld 
regime of lubrication for the general case of elllptlcalcontacts. In this 
regime several formulas of the lubricant film thickness have been proposed by 
Hamrock and Dowson, by Dowson et al., and more recently by Houpert. However, 
either they do not Include the load parameter W, which has a strong effect on 
film thickness, or they overestimate the film thickness by using the Barus 
formula for pressure-viscosity characteristics. In the current study the 
Roelands formula has been used for the pressure-viscosity relationship. The 
effects of the dimensionless load, speed, and materials parameters, the radius 
ratio, and the lubricant entrainment direction have been Investigated. The 
dimensionless load parameter was varied over a range of one order of magnitude 
The dimensionless speed parameter was varied by 5.6 times the lowest value. 
Conditions corresponding to the use of solid materials of steel, bronze, and 
Silicon nitride and lubricants of paraffinic and naphthenic mineral oil were 
considered In obtaining the exponent In the dimensionless materials parameter. 
The radius ratio was varied from 0.2 to 64 (a configuration approaching a line 
contact). Forty-one cases were used In obtaining the minimum film thickness 
formula: 

H - 1 78G^’ 386yl • 266^-0.880^ _ g -0.0387aj 
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Contour plots Indicate In detail the pressure developed between the contacting 
solids . 

NOMENCLATURE 

D difference, [(Hq - Hq)/Hq]X 100, percent 

2 

E modulus of elasticity, N/m 

E' effective elastic modulus, 2/[(l - v 2 )/E A + (1 - » 2 )/Eg], N/m 2 

G dimensionless materials parameter, E'/p. 

K 1v,as 

H dimensionless film thickness, h/R 

H , dimensionless minimum film thickness, h , /R 

min min x 

H q dimensionless central film thickness, h^/R 

Hp calculated dimensionless central film thickness from least-squares 

analysis 

h f 1 lm thickness , m 

h , minimum film thickness, m 

min ' 

h Q central film thickness {minimum film thickness as well In 

pi ezovl scous-rigld lubrication regime), m 

P dimensionless pressure, p/E * 

p pressure. Pa 

2 

p, asymptotic Isovlscous pressure, N/m 

I V j(l j 

R effective radius, m 

r radius ratio 

s,t constants defining fluid, used In Eq. (8) 

U dimensionless speed parameter, ti u/E'R 

n — ~ x 

u yu + u , m/s 

T x y ’ 

surface velocity In x direction, m/s 

Uy surface velocity in y direction, m/s 

2 

W dimensionless load parameter, w/E'R x 
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w 

X.Y 

x 

y 

z 

Ct 

n 

n 

n 0 

e 

V 


normal applied load, N 

dimensionless coordinates, x/R , y/R 

x x 

coordinate In rolling direction, m 

coordinate In transverse direction, m 

viscosity pressure Index, a dimensionless constant 

radius ratio, R /R 
y x 

lubricant viscosity, (N sj/m? 
dimensionless lubricant viscosity, n/ng 

2 

lubricant viscosity at atmospheric pressure, (N s)/m 

angle between lubricant entrainment vector and rolling direction, deg 

Poisson's ratio 


p lubricant density, kg/rn^ 

p dimensionless density, p/pg 

3 

p 0 lubricant density at atmospheric pressure, kg/m 

3/2 

<p PH 

Subscripts 
a solid a 

b solid b 

x,y coordinates In plane of lubricating film 

INTRODUCTION 

The development of fluid film lubrication, Implying complete separation 
of the surfaces and no asperity contact, Is Influenced by two major physical 
effects: the elastic deformation of the solid under an applied load, and the 

variation of fluid viscosity with pressure. Therefore It Is possible to have 
four main regimes of fluid film lubrication (Hamrock and Dowson, 1981) 
depending on the magnitude of these effects and on their Importance. The 
regimes can be described as (1) 1 sovl scous-rlgld , (2) plezovl scous-rlgid , 

(3) Isovlscous-elastic, and (4) plezovl scous-elastlc . 
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The hydrodynamic lubrication (Isovlscous-rlgld) theory was first applied 
to nonconformal contacts by Martin (1916) to explain the mechanism of gear 
lubrication. His solution of the Reynolds equation for the lubrication between 
two rigid circular cylinders by means of an Isovlscous, Incompressible fluid 

t 

presented a useful beginning to theoretical studies. Nevertheless Martin's 
theoretical work discouraged the view that spur gears could be lubricated by 
hydrodynamic action since his results Indicated that the film thickness of 
spur gear lubrication was quite small compared with the surface roughness. A 
remarkable extension of the classical analytical solution for elliptical 
contacts lubricated by either an Isovlscous or plezovlscous fluid was obtained 
by Kapltza (1955). However, applying the half-Sorranerfeld boundary condition 
used In Kapltza's analysis violated flow continuity at the cavitation boundary. 

Brewe et al. (1979) obtained a film thickness equation for the lubrication 
of fully flooded, rigid, Isovlscous point contacts through a numerical . analysis 
that used a more realistic Reynolds boundary condition for the film rupture In 
the exit region. He found that the minimum film thickness had the same speed, 
viscosity, and load dependence as In Kapltza's classical solution. However, 
the Incorporation of the Reynolds boundary condition resulted in an additional 
geometry effect. 

In the plezovlscous-rlgld regime of lubrication the pressure within the 
contact Is sufficiently high to Increase the fluid viscosity significantly 
within the contact while the deformation of the surfaces remains an 
Insignificant part of the fluid film thickness. This form of lubrication Is 
encountered on roller end-guide flanges. In contacts In moderately loaded 
cylindrical tapered rollers, and between some piston rings and cylinder liners. 

In the plezovlscous-rlgld regime of lubrication Marko and Clegg (1979) 
proposed a formula obtained by means of curve fitting Dowson and Whitaker's 
(1965) results for line contacts. For point contacts Hamrock and Dowson 
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(1978) proposed as an Interim measure that Blok's (1952) solution for line 
contacts be adjusted by applying the same "side-leakage factor" as that 
derived for the plezovlscous-elastlc condition. Oowson et al. (1983) obtained 
the numerical solution for a compressible Newtonian lubricant exhibiting 
pressure-viscosity characteristics and subject to the Reynolds cavitation 
boundary condition. The formulas from the last two papers describe the 
"limiting" film thickness generating Infinite pressures, as discussed by Blok 
(1952). However, they do not Include the load parameter W, which has a 
strong effect on film thickness when plezovlscous effects are considered. 
Houpert (1984) developed a sophisticated general formula by curve fitting 
Dalmaz's (1979) results. A shortcoming of Dalmaz's (1979) work Is that he 
used the Barus exponential formula for pressure-viscosity characteristics. 

That formula tends to give higher values of viscosity than the results 
obtained from Roelands formula (Jones et al., 1975). Furthermore Dalmaz's 
(1979) results for the Isovlscous case produced a lower exponent on W/U than 
did Brewe et al . (1979). This appears to be due to starvation effects entering 
the fully flooded results from the designation of the Inlet boundary condition. 

Figure 1, obtained from Meuleman et al. (1985), shows the differences 
between the Barus formula, the Roelands formula, and the experimental data of 
Hirst and Moore (1979). It Is apparent that the Roelands formula represents 
the experimental result for the pressure-viscosity relationship more accurately 
than the Barus formula at high pressures. Therefore the Roelands formula Is 
used In this study. The researchers In the past used the Barus formula and 
neglected the lubricant compressibility to obtain straightforward numerical 
analysis of the resulting linear partial differential equation. The numerical 
analysis employed In this study Is more complicated because the resulting 
partial differential equation Is now nonlinear. 
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In the current study the numerical solution for the plezovlscous-rlgld 
lubrication regime Is presented for the full spectrum of conditions. The 
Influence of lubricant entrainment direction has also been studied. The 
condition that the lubricant entraining vector Is not parallel to the minor 
axis of the contact might arise, for example. In helical, spiral bevel, and 
hypold gears. The radius ratio was varied from 0.2 to 64 to cover any contact 
ranging from something similar to a disk rolling on a plate (a < 1) to a 
contact approaching a nominal line contact (a -» 70) such as a barrel-shaped 
roller against a plate. The effects of the dimensionless speed, load, and 
materials parameters on minimum film thickness were also Investigated. The 
dimensionless load parameter was varied over a range of an order of magnitude. 
The dimensionless speed parameter was varied over a range 5.6 times the lowest 
speed value. Conditions corresponding to the use of solid materials of steel, 
bronze, and silicon nitride and lubricants of paraffinic and naphthenic mineral 
oils were considered In obtaining the exponent In the dimensionless materials 
parameter. Forty-one cases were used to obtain a simple empirical minimum 
film thickness formula. Contour plots are shown that Indicate In detail the 
pressure developed between the solids. 

METHOD OF CALCULATION 
Reynolds Equation 

For the coordinate systems x,y the general Reynolds equation for point 
contact can be expressed as 


a 

eh! a£| 
Vl2n 3 x) 

. a_ 

Vl2n ayy 

a_ 

>(u + u. )h" 

ax bx 7 

. a_ 

p ,(u av * VI 

ax 

ay 

ax 

2 

ay 

2 


where 

u surface velocity of solid a In x direction 

d X 

u bx surface velocity of solid b In x direction 
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u ay surface velocity of solid a In y direction 

u^ surface velocity of solid b In y direction 

If the surface velocities are assumed to be constant, the Reynolds equation 
can be expressed as 



u iL(ph) 

x ax 


+ u 


3_(pb) 

y ay 


(2) 


where 
u = 


(u ♦ u K )/2 
ax bx 7 


u = (u + u. )/2 
y ' ay by 7 

By Introducing U and e, where 
Eq. 2 becomes 


u = yu* . u 


and e = tan 


-1 


(u /u ) 
' x y 7 


_ 12u [cos e - a lp- h > + sin 0 

3x \ n 3x / + 3y\ n 3y/- 12u [ cose 3x + s1n 6 3y j 


a(ph) I 


(3) 


By letting 



Eq. (3) can be described In the dimensionless form as 
3 3P \ A 3 


fy) * 1?u cos 6 m 


+ sin e 


iLpHlI 
3Y J 


(4) 


where 

u = n n u/E 1 R 
0 x 

Equation (4) Is the Reynolds equation for which the dimensionless pressure P 
will be determined. However, before proceeding, the dimensionless density 
pg, the dimensionless viscosity ng, and the dimensionless film thickness 
H will have to be expressed. 
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Pressure-Viscosity Formula 

It Is generally known that the viscosity of a lubricant Is a function of 

• -* t f • 

temperature and pressure. A generally accepted relationship that describes 
the pressure-viscosity dependency Is the Barus equation 



where 

n viscosity at gauge pressure 

tiq viscosity at atmospheric pressure 

a pressure-viscosity coefficient of lubricant 

Unfortunately pressure-viscosity data seldom follow this simple relationship, 
and It Is valid as a reasonable approximation only In a moderate-pressure 
range. 

Roelands (1966) In a more extensive experimental study of the effect of 
pressure on the viscosity of lubricants has developed an empirical formula 
written as 

log n + 1.200 = (log n 0 + 1 - 200) ^1 + 

where 

p gauge pressure, kgf/cm 2 

Z viscosity Index, a dimensionless constant 

Rearranging terms gives 

5 ' (l+p/2000) 2 • in 1.2[( Up/2000) Z -l] (7) 

n = t)q x i u 

The temperature effect Is normally accounted for In Hq. 

Pressure-Density Formula 

The variation of density with pressure Is roughly linear at low pressures, 
but the rate of Increase falls away at high pressures. From Dowson and 
Hlgglnson (1966) the dimensionless density for mineral oil can be written as 
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= 1 + 


J_E_ 


( 8 ) 


P 0 1 + tp 

where s and t are constants that depend on the fluid. 

Film Thickness 

The separation of two rigid solids a and b having radii of curvature 
r ,r and r. .r. * he vicinity of the point of closest approach can be 

aX aj DX Dy 

considered as a geometrically equivalent solid of principal radii R ,R 

x y / 

adjacent to a plane, as shown In Fig. 2. The effective radius of curvature can 
be expressed as 


1 1 

r r. 
ax bx 


1 ^ 1 
r ay r by 


(9) 


It Is assumed that convex surfaces exhibit positive curvature, and concave 
surfaces negative curvature. Therefore, If the center of curvature lies within 
the solid, the radius of curvature Is positive; If the center of curvature lies 
outside the solid, the radius Is negative. 

The film thickness between two rigid bodies In point contact can be 
written as 

h = h Q * S(x,y) (10) 

The separation of two rigid solids (Fig. 2(a)) In which the principal axes of 
Inertia of the two bodies are parallel can be written as 


S = S + S. + S + S. 
ax bx ay by 
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where 


= r - V r 
ax t 

-{ 

-V 

by " i 


ax 


$bx ~ r bx 


S = r 
ay ay 


2 

ax 


r bx _ x 


2 2 
r ay " y 


2 2 
- y 


(11) 


by " by ▼ by 

The separation In terms of the coordinates and the effective radius of 
curvature (Fig. 2(b)) Is 


S(x,y) = R x - 


2 2 

R x ' x + R y - 




D 2 2 

Ry - y 


( 12 ) 


Substituting Eq. (12) Into Eq. (11) while at the same time making this 
equation dimensionless gives 

r t — 

,2 


H = H, 


♦ ' - Vi - X 2 -(;/ 


(13) 


where a = R /R . 

y x 

Boundary Conditions 

There are two boundary conditions: 

(1) The Reynolds boundary condition Is applied In the divergent film. 
Namely, P = dP/dX = dP/dY = 0 at the cavitation boundary. 

(2) The pressure on the boundaries of the computation region Is assumed 
to be zero. The fully flooded conjunction Is considered In the present study 
( 1 - e . , Increasing the computation region of the conjunction does not affect 
the results). 

NUMERICAL ANALYSIS 


Having defined the density, viscosity, and film thickness, we are In a 
position to solve the Reynolds equation subjected to the Reynolds cavitation 
boundary condition. The dimensionless pressure P Is notorious for Its steep 
gradient, which Is not welcome when performing numerical analysis via 
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relaxation methods. To produce a more gentle curve, a parameter <p Is 
Introduced where 

3/2 

<p = PH J,C (14) 

This substitution also has the advantage of eliminating all terms containing 
derivatives of products of H and P or H and <p. Substituting Eq. (14) 
Into Eq. (4) , we obtain 


j3/2 




1/2 3H 
ax 


+ H 


3/2 



3 3 u l/2 3H 

2 » av ( ? H Sy 


= 12 cos e 


a(pH) 

ax 


+ sin e 


a(p.H) ~ 
aY 


(15) 


The method of frozen coefficients (sometimes also referred to as Kacanov's 
method (Fuclk et al . (1975)) Is applied to solve this nonlinear partial 
differential equation. Second-order central finite difference approximations 
to the modified form of the Reynolds equation are used, thus forming a set of 
algebraic equations that are solved by the Gauss-Seldel Iterative method with 
overrelaxatlon. For optimal efficiency a variable-mesh structure Is used to 
enhance the accuracy In the region of high pressure and large pressure 
gradient. 

The normal applied load can be evaluated by 

T y m f\ 

w = I I p(x,y) dx dy (16) 

J *o Jx o 


The double Simpson's Integration technique Is applied to the Integration. 


RESULTS 


Dimensionless Grouping 

From the variables of the numerical analysis the following dimensionless 

groupings are written: 

(1) Dimensionless film thickness 

H=h/R x 
11 


( 17 ) 



(2) Dimensionless load parameter' 

W = w/E'R 2 (18) 

(3) Dimensionless speed parameter 

U = n n u/E , R y (19) 

V 2 . 2 • 

u "+‘-u . • ’■ 

^ j 

(4) Dimensionless materials parameter 

6 = £,/ Va S (20 > 

where p^ v Is the asymptotic, Isovlscous pressure gradient obtained from 
Roelands (1966) 

(5) Radius ratio 

a = R y /R x (21) 

(6) Lubricant entraining angle 

e = tan _1 (u /u ) In degrees (22) 

y * 

The dimensionless film thickness can thus be written as a function of the 
other five parameters: 

Hg = f(W,U, 0,6,6) (23) 

The most Important practical aspects of hydrodynamic lubrication of 

nonconformal contacts Is the determination of the minimum film thickness 

within the contact. Therefore In the fully flooded results to be presented 

the dimensionless parameters (W, U, o, G, and 6) will be varied, and the 

effect on the minimum film thickness will be studied. Note that In Eqs. (17) 

to (21), by changing the normal applied load w, the dimensionless load W Is 

changed, but the other dimensionless parameters remain unchanged. Similar 

statements can be made about changing u^ In the dimensionless speed U and 

R In the radius ratio, 
x 

Influence of Load 

The dimensionless parameters U, G, a, and e were kept constant while 
varying the dimensionless minimum film thickness Hg to get the dimensionless 
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load capacity W at each different Hq. The values at which the remaining 
parameters U, 6, a, and e were held constant during the calculations 
were U = 0.16833xl(f 11 . G = 4522, a = 16, and 6 = 0. 

Table 1 shows the computed load capacities for 10 values of minimum film 
thickness. From these 10 pairs of data It Is possible to determine a good 
empirical relationship between the minimum film thickness Hq and the load 
capacity W: 

c 

«0 = C-|W (24) 


By applying a least-squares power fit to the 10 pairs of data [(W. , 

Hq i ) , 1 = 1,2 10], the values of C 1 and C ? were found to be 

C.J = 2.60615x10 12 « 2.606xl0 -12 and = -0.88019 « -0.880. From the value 
of C 2 and Eq. (24) the effect of load on minimum film thickness Is written as 


H„ ■« W 


■0.880 


(25) 


In Fig. 3 the variation of dimensionless minimum film thickness with 
dimensionless load Is plotted for 10 data points. On a log-log scale the plot 
Is linear. The percentage difference between the minimum film thickness 
obtained from computational results Hq and the minimum film thickness 
obtained from the least-squares fit equation Hq Is expressed as 

0 “ [<»o " V /H 0> X10 ° 

Influence of Speed 

If the surface velocity u Is changed, the dimensionless speed parameter 
U Is modified as shown In Eq. (19), but the other dimensionless parameters 
(Hq, G, a, and 6) remain constant. The values at which these dimensionless 
parameters were held constant In the calculations performed to determine the 
Influence of speed on film thickness are 

a = 16, Hq = 4.8xl0“ 6 , G = 4522, 6=0 (27) 
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Values of the dimensionless speed parameter U and the corresponding 
dimensionless load capacities as obtained from the numerical computations are 
presented In Table 2. Since the relationship between the minimum film 
thickness Hg and the load capacity W has been obtained, the relationship 
between the minimum film thickness and the speed parameter can be written In 
the form 


' H 0 /W - 0 - 880 



(27) 


By applying a least-squares fit to the six pairs of data [U^, Hg , 

1 = 1 , . . . , 6 ] , the values of C 3 and were found to be 

C 3 = 1.4899X10" 13 » 1.49X10" 13 and C 4 = 1.2655 a 1.266. From the 

value obtained for C. and Eq. (27) the effect of dimensionless speed on 

4 

dimensionless film thickness can be written an 

H 0 “ y1 266 (28) 
The data presented In Table 2 are plotted In Fig. 4. As was true for the 

load-versus-f 1 lm-thlckness results, the speed-versus-f 1 Im-thlckness results 

are linear when plotted on a log-log scale. 

Influence of Radius Ratio 

To study the effect of geometry on minimum film thickness, the 
dimensionless parameters for minimum film thickness Hg, speed U, materials 
G, and lubricant entrainment direction e were held constant at the 
following values: 

Hg = 4.8xl0 -6 , U = 0.1 683x10” 11 , G = 4522, 0=0 
Table 3 shows the computed values of dimensionless load capacity for 22 
values of radius ratio. It Is possible to determine a good empirical 
relationship between Hg and a for the condition considered In the 
computation. The form of relationship chosen after Investigating a number of 
different expressions can be written as 
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1 - (H 0 /H 0r ) - Ae B 


(29) 


... -0.0387a. 

H q * (1 - e ) 


where H n Is chosen to be the film thickness at rectangular contact. 

A least-squares exponential curve was fitted to the 22 pairs of data 
points to obtain values for A and B In Eq. (29): A = 0.989 * 1.00 and 
B = 0.03866 « 0.0387. Substituting these values Into Eq. (29) resulted In 
the following relationship between the radius ratio and minimum film thickness: 

(30) 

The effect of radius ratio on film thickness for the 22 data points Is shown In 

Fig. 5. It Is most significant that the computed value of A Is approximately 

unity since this ensures that the minimum film thickness approaches zero as 

the radius ratio goes to zero. Contour plots for three radius ratios (l.e., a 

of 16, 8, and 4) are shown In Fig. 6 to demonstrate the Influence of geometry. 

Since the Isobars In each case are evenly spaced, the pressure gradients can 

be easily depicted. As the radius ratio Increases, the steeper pressure 

gradients are predominantly along the rolling direction. This Implies that 

the amount of side leakage decreases as radius ratio Increases. Figure 7 shows 

that the geometry effects In four regimes of lubrication (l.e., hard EHL 

(Hamrock and Dowson, 1977), soft EHL (Hamrock and Dowson, 1978), 

plezovlscous-rlgld, and Isovlscous-rlgld (Brewe et al., 1979)), where H m ^ n r 

Is the minimum film thickness for rectangular contacts. It Is assumed that 

when a = 150, the rectangular contact limit Is realized. The ratio 

H . /H . approaches the limiting value very quickly In hard EHL, and 
min min, r 

most slowly In the Isovlscous-rlgld case. 

Influence of Material Properties 

A study of the Influence of the dimensionless materials parameter 6 on 
minimum film thickness has to be approached with caution since In practice It 
Is not possible to change the physical properties of the materials, and hence 
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the value of 6, without influencing the other dimensionless parameters 
considered earlier. Equations (18) to (20) show that when either the materials 
of the solids (as expressed In E 1 ) or the other lubricants (as expressed In 
n» and p. ) are varied, not only does 6 change, but so do the 
dimensionless speed U and load W parameters. Only the radius ratio, 
minimum film thickness, and lubricant entrainment direction can be held fixed; 
and for all results presented In Table 4, Hq = 4.8xl0~ & or 1.4xl0~ 6 , a = 16, 
and e = 0 are adopted. 

The results obtained from calculations performed for six values of the 
dimensionless materials parameter are summarized In Table 4. The general form 
of these results, showing how the minimum film thickness Is a function of the 
dimensionless materials parameter, Is 

c ft 

C = C ? G (31) 

where 


C = H q /[(1 


-0.0387a. „1 .266^,-0.880. 
e ) U W ] 


By applying a least-squares power fit to the six pairs of data, the values of 
C and C Q were found to be C = 175.04 and C = 0.386096 « 0.386. 

The effect of the dimensionless materials parameter on the minimum film 
thickness Is approximated with adequate accuracy as 


,0.386 


(32) 


The variation of dimensionless materials parameter G with dimensionless film 
thickness H Q Is shown In Fig. 8. The six data points given In Table 4 are 
shown In this figure. Again, on the log-log scale the representation of the 
materials parameter versus film thickness Is linear. 

Influence of Lubricant Entrainment Direction 

If the velocity of the entrainment lubricant Is kept constant but the_ 
component velocity In the x direction u^ and the component velocity In the 
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y direction u y are changed, the lubricant entraining angle Is modified as 
shown In Eq. (22) while the other dimensionless parameters (Hq, G, and a) 
remain constant. The values at which these dimensionless parameters are held 
constant In the calculations performed to determine the effect of lubricant 
entrainment direction are 

a = 20, H q = 5.6x1 0 -6 , G = 4522, U = 0.75749X10 -11 
The results obtained from calculations performed for 10 values of the 
lubrication entrainment direction are summarized In Table 5. A cosine function 
was chosen for curve fitting, namely 

T = C g e (33) 

where T = cos -1 (W”^'® 80 /Wq®‘® 8 ^) and Wq denotes the load capacity at 
e=0. By applying a least-squares fit to the 10 pairs of data, the value 
of Cg was found to be C g = 0.0200027 « 0.0200. Therefore the effect 
of the lubricant entrainment direction on the minimum film thickness can be 
written as 

H q <* cos (1 .1460) (34) 

This Is shown In Fig. 9. It Is significant that when the entraining angle Is 
zero, the cosine function turns out to be 1 . In other words, there Is no 
effect of lubricant entrainment direction on minimum film thickness. Although 
the Influence of lubricant entrainment direction was Investigated for 0 = 20° 
only. It can be applied to the cases from © = 16° to 0 = 30°. 

Minimum Film Thickness Formula 

The proportionality Eqs. (25), (28), (30), and (34) have established how 
the minimum film thickness varied with the dimensionless load, speed, and 
materials parameters and radius ratio, respectively. This enables a composite 
dimensionless minimum film thickness formula for a fully flooded, rigid, 
elliptic-contact solid lubricated by a plezovlscous fluid to be modeled by 
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(35) 


u , __ _0 . 386. .1 .266,. -0.880/, -0.0387a \ 

Hq = i /ob u w - e / 

with 85.6 percent confidence (except for the small value of a) ±8.92 percent 
of the value from the numerical analysis. Equation (35) can be written In 
dimensional form as 


h na/i/n xO.386. J.266 u -0.880 d 1.494 

h 0 " 178(1/p 1v,as ) (Un 0 ) W R 




-0.0387R /R > 
y x 


(36) 


Note that the effective elastic modulus E' does not appear In this equation. 
Cases 1 to 41 In Table 6 were used In obtaining Eq. (35). To gain more 
confidence In the application of this empirical equation, some other values of 
dimensionless speed U and minimum film thickness Hq not In the range for 
curve fitting were checked. The results are shown In cases 52 and 53 of Table 
6. The good agreement between the prediction from the minimum film thickness 
formula and computed results Is encouraging. 


DISCUSSION 

A numerical solution for the plezovlscous effect In nonconformal rigid 
contacts lubricated hydrodynamlcally has been developed. The following remarks 
can be made: 

(1) The effective elastic modulus E' does not appear In the dimensional 
film thickness equation. This corresponds to a rigid contact; there Is no 
effect of elastic properties. 

(2) In all cases. If the maximum dimensionless viscosity was less than 3, 
the load capacity was no more than 12 percent larger than the load capacity 
without the plezovlscous effect. This Indicates a very small plezovlscous 
effect. All the data sets used to generate Eq. (35) had maximum dimensionless 
viscosity larger than 3. 
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(3) The exponents of dimensionless load and dimensionless speed In the 
minimum film thickness equation for the plezovlscous-rlgld regime were between 
those for the plezovlscous-elastlc regime and the Isovlscous-rlgld regime. 

(4) The film thickness formula according to the side-leakage factor 
proposed by Hamrock and Dowson (1978) Is 

H q = 1.66(GU) 2/3 {1 - exp[-0. 68(1 .03a 0,66 )]} (37) 

The more recent formula from the numerical solution by Dowson et al. (1983) Is 

H 0 = 1.66(GU) 2/3 ( T -^|' 6 (38) 

Figure 10 compares the minimum film thickness as obtained from Eqs. (37) and 
(38) with the present result for W = 0.6x10 7 or 0.9x10 7 and a = 8 or 16. 

The deviations are large for small values of GU. The deviations resulted 
because the load parameter W was not Included In Eqs. (28) and (29), giving 
erroneous results when the plezovlscous effects are small (l.e., GU Is small), 
and because the Barus formula was used for the pressure-viscosity correlation 
In Eqs. (37) and (38) and caused an overestimate of the minimum film thickness. 

(5) The film thickness formula proposed by Houpert (1984) by means of 
curve fitting the numerical solution of Dalmaz Is 



with the formula proposed by Houpert. Table 7 shows the ratio of 

(Eq. (39)) to Hq (Eq. (35)). The deviation Is between 28 and 54 percent. 
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A possible cause of this difference may be that Houpert's (1984) results Imply 
an Incompressible fluid and the Barus pressure-viscosity formula. 

SUMMARY OF RESULTS 

A procedure for the numerical solution of the complete plezovlscous-rlgld 

regime of lubrication for nonconformal contacts has been demonstrated. In the 

3/2 

numerical analysis of the Reynolds equation, the parameter <p = PH was 
Introduced to help the relaxation process. A variable-mesh nodal structure 
was used to provide close spacing In and around the pressure peak. A more 
realistic formula Is used for the pressure-viscosity relationship. Lubricant 
compressibility was also considered. 

By using the procedures outlined In the analysis, the effects of the 
dimensionless load W, speed U, and materials G parameters, radius ratio 
a, and lubricant entraining direction e on minimum film thickness were 
Investigated. Forty-one cases were used to generate the minimum film thickness 
relationship 

u ., Qr 0.386 M l .266. .-0.880.. -0. 0387a. 

Hq = 1 78G U W ( 1 - e ) 

The most dominant exponent occurs In association with the speed 
parameter; the exponent on the load parameter -0.880 Is between -2 for the 
Isovlscous-rlgld regime and -0.073 for the plezovlscous-elastlc regime. The 
materials parameter also carries a significant exponent, although the range of 
the parameter In engineering Is limited. The geometry effect shows the same 
curve as In elastohydrodynamlc lubrication but approaches the limiting value 
much more slowly. The effect of the direction of lubricant entrainment Is a 
cosine function for the geometries studied. The dimensionless load parameter 
values used In obtaining the preceding equation were varied over a range of 
one order of magnitude. The dimensionless speed values were varied by 5.6 
times the lowest value. Conditions corresponding to the use of solid 
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materials of steel, bronze, and silicon nitride and lubricants of paraffinic 
and naphthenic mineral oils are considered In obtaining the exponent In the 
dimensionless materials parameter. The radius ratio was varied from 0.2 to 64 
(a configuration approaching a line contact). Contour plots are shown that 
Indicate In detail the pressure developed between the solids. 

REFERENCES 

Blok, H., "Discussion of Paper by E. McEwen," 3. Inst. Petrol. . 38, 673-683, 
(1952). 

Brewe, D.E., Hamrock, B.J., and Taylor, C.M., "Effect of Geometry on 
Hydrodynamic Film Thickness," J. Lubr. Techno!, . 101 , 231-239, (1979). 
Dalmaz, G., "Le film mince vlsqueux dans les contacts hertzlens en regimes 
hydrodynamlque et elastohydrodynamlque," Docteur d'Etat Es Sciences Thesis, 
I.N.S.A. Lyon, (1979). 

Dowson, D. and Hlgglnson, G.R., Elastohydrodynamlc Lubrication . Pergamon 
Press, (1966). 

Dowson, D. and Whitaker, A.V., "The Isothermal Lubrication of Cylinders," ASLE 
Trans. . 8, 224-234, (1965). 

Dowson, D., Dunn, J.F., and Taylor, C.M., "The Piezo-viscous Fluid, Rigid 
Solid Regime of Lubrication," Proc. Inst. Mech. Eng. , 197c , 43-52, (1983). 
Fuclk, S., Kratochvll, A., and Necas, 3., "Kacanov's Method and Its 
Application," Rev. Roum. Hath. Pures Appl. . 20, 907-916, (1975). 

Hamrock, B.3. and Dowson, D., "Isothermal Elastohydrodynamlc Lubrication of 
Point Contacts, Part III - Fully Flooded Results," 3. Lubr. Technol . . 99, 
264-276, (1977). 

Hamrock, B.3. and Dowson, D., "Elastohydrodynamlc Lubrication of Elliptical 
Contacts for Materials of Low Elastic Modulus, Part I - Fully Flooded 
Conjunction," J. Lubr. Technol. . 100 . 236-245, (1978). 


21 


Hamrock, B.J. and Dowson, 0., Ball Bearing Lubrication - The 

Elastohydrodynamlcs of Elliptical Contacts . John Wiley & Sons, New York, 
(1981). 

Hirst, W. and Moore, A.J., "Elastohydrodynamlc Lubrication at High 

Pressure, Part 2 - Non-Newtonian Behavior," Proc. R. Soc.. London. A. . 365 . 
537-565, (1979). 

Houpert, L., "The Film Thickness In Plezovlscous-Rlgld Regime; Film Thickness 
Lubrication Regimes Transition Criteria," J. Lubr. Technol. . 106 , 375-385, 
(1984). 

Jones, W.R., Jr., Johnson, R.L., Winer, W.O., and Sanborn, D.M., 

g 

"Pressure-Viscosity Measurements for Several Lubricants to 5.5x10 Newtons 
per Square Meter (8xl0 4 psl ) and 149 °C (300 °F)," ASLE Trans. . 18. 

249-262, (1975). 

Kapltza, P.L., "Hydrodynamic Theory of Lubrication during Rolling," Zh. Tekh. 
Flz. . 25, 747-762, (1955). 

Markho, P.H., and Clegg, D.B., "Reflections on Some Aspects of Lubrication of 
Concentrated Line Contacts," J. Lubr. Technol. . 101 . 528-531, (1979). 

Martin, H.M., "The Lubrication of Gear Teeth," Engineering (London) . 102 . 
119-121, (1916). 

Meuleman, P.K., Lubrecht, A. A., and Napel, W.E., "Traction In 

Elastohydrodynamlc Lubrication - High Contact Pressure Experiments and a 
Thermal Roelands-Eyrlng-Maxwell Model," Department of Mechanical Engineering, 
Twente University of Technology, Enschede, The Netherlands (1985). 

Roelands, C.J.A., Correlational Aspect of the Viscosity-Temperature-Pressure 
Relationship of Lubricating Oils ." Druk V.R.B., Grlnlngeon, Netherlands, 
(1966). 


22 



TABLE 1 - EFFECT OF DIMENSIONLESS LOAD PARAMETER ON 
MINIMUM FILM THICKNESS 


Dimensionless 
load parameter, 
U 

Minimum film thickness 

Difference, 

D, 

percent 

Obtained from 
computational 
result, 

H 0 

Obtained from 
least-squares 
fit, 

H q 

0 .90971xl0“ 7 

4.1xl0“ 6 

4 .13xlO -6 

0.77 

.89091 

4.2 

4.21 

.20 

.85159 

4.4 

4.38 

-.48 

.80860 

4.6 

4.78 

-.37 

. 77376 

4.8 

4.76 

-.75 

.74467 

5.0 

. 4.93 

-1.45 

.68766 

5.2 

5.29 

1.64 

.65859 

5.5 

5.49 

-.18 

.64956 

5.6 

5.56 

-.77 

.63293 

5.8 

5.69 

-1.98 . 


TABLE 2 - EFFECT OF DIMENSIONLESS SPEED PARAMETER ON MINIMUM FILM THICKNESS 


Dimensionless 
speed parameter, 
U 

Dimensionless 
load parameter, 
U 

Minimum film thickness 

Difference, 

D,. 

percent 

Obtained from 
computational 
result, 

H 0 

Obtained from 
least-squares 
fU, 

H 0 

0.21883xl0 -11 

SI 

4.8 

4.68 

-2.56 

.20120 

• ' • 1 



4.83 

■ .72 

.18516 

.88950 



4.84 

.76 

.16833 

.77377 



4.85 

.96 

.15150 

.67629 



4.77 

-.53 

.13466 

.58454 



4.68 

-2.58 


















TABLE 3 - EFFECT OF ELL IPTIC ITY PARAMETER ON MINIMUM FILM 

THICKNESS 



TABLE 4 - EFFECT OF DIMENSIONLESS MATERIALS PARAMETER ON MINIMUM FILM THICKNESS 


Solid 

material 

Lubricant 

Dimensionless parameter 

Minimum film thickness 

Difference, 


Materials, 

G 

Speed, 

U 

Load, 

W 

Obtained from 
computational 
result, 

H 0 

Obtained from 
least-squares 
fit, 

H 0 

D, 

percent 

Steel 

Paraffinic 

4 522 

0.16832xl0” u 

0.77376x10”^ 

4.8xl0” 6 

4.68xl0” 6 

-2.48 

Bronze 

Paraffinic 

2 310 

.32957 

1.51499 

4.8 

4.682 

-2.45 

Si 1 icon 
nitride 

Paraffinic 

6 785 

.11218 

.51571 

4.8 

4.682 

-2.45 

Steel 

Naphthenic 

7 031 

.48123 

1.18602 

1.4xl0” 5 

1.441 

2.92 

Bronze 

Naphthenic 

3 591 

.94227 

2.32216 

1.4 • 

1.441 

2.92 

Si licon 
nitride 

Naphthenic 

10 549 

.32073 

.79047 

1.4 

1.441 

2.92 



























TABLE 5 - EFFECT OF LUBRICANT ENTRAINMENT DIRECTION ON MINIMUM FILM 

THICKNESS 


Luoricant . 
entrainment 
direction , 
9, 
deg 

Dimensionless 
load parameter, 
U 

Minimum film thickness 

Difference 

0, 

percent 

Obtained from 
computational 
result, 

H 0 

Obtained from 
least-squares 
fit, 

H 0 

0 

0.l4965xl0“ 6 

5.6xl0“ 6 

5.60x10“® 

0 

5 

.14841 



5.56 

-.69 

10 

.14532 



5.54 

-1.03 

15 

.14043 



5.52 

-1.49 

20 

.13427 



5.50 

-1.80 

25 

.12708 



5.50 

-1.80 

30 

.11941 



5.54 

-1.14 

35 

.11902 



5.60 

0 

40 

.10200 



5.49 

-1.94 

45 

.09278 

’ 


5.33 

-4.88 











TABLE 6 - DATA SHOWING EFFECT OF LOAD, SPEED,' RADIUS RATIO, MATERIALS PARAMETER, AND LUBRICANT ENTRAINMENT DIRECTION ON 

MINIMUM FILM THICKNESS 


Case 

Dimensionless 

load 

parameter, 

U 

Dimensionless 

speed 

parameter, 

U 

Radius 

ratio, 

a 

Dimensionless 

materials 

parameter, 

G 

Lubricant 

entrainment 

direction, 

8, 

deg 

Minimum film thickness 

Difference, 

D, 

percent 

Results 

Obtained from 
computational 
results, 

H 0 

Obtained from 
least-squares 
fit, 

H 0 


0.90971x10“^ 

0.16833xl0“ U 

16 

4 

522 

0 

4 

IxlO - 6 

4.13xl0“ 6 

0.73 

Load 


.89091 









4 

2 

4.21 

.24 




.85159 









4 

4 

4.38 

-.47 




.80860 









4 

6 

4.58 

-.36 




.77376 









4 

8 

4.76 

-.74 




.74467 









5 

0 

4.93 

-1.44 




.68766 









5 

2 

5.11 

-1.65 



8 

.65859 









5 

5 

5.49 

-.17 



9 

.64956 









5 

6 

5.56 

-.76 



10 

.63293 









5 

8 

5.69 

1.97 



11 

1.17505 

.21883 







4 

8 

4.60 

-4.20 

Speed plus 

12 

1.00851 

.20120 









4.75 

-.97 

case 5 

13 

.88950 

.18516 









4.35 

-.94 



14 

.67629 

.15150 









4.69 

-2.20 



15 

.58454 

.13466 









4.60 

-4.22 



16 

1.51499 

.32957 



. 2 

310 





4.64 

-3.33 

Materials 

17 

.51571 

.11218 



' 6 

785 





4.64 

-3.33 

plus case 5 

18 

1.18602 

.48123 



7 031 



14xl0“ b 

1 .38x10“ 

-2.14 



19 

2.32216 

.94227 



3 

591 



14 


1.38 

-2.14 



20 

.79047 

.32073 



10 

549 



14 


1.38 

-2.14 



21 

1.61064 

.16833 • 

64 

4 522 



4.8xl0“ b 

4.96xl0“ 6 

3.32 

Radius 

rati 

22 

1.59341 



61 







4.95 

3.13 

plus 

case 

23 

1.56632 



59 







4.98 

3.82 



24 

1.55278 



58 







4.99 

4.14 



25 

1.49756 



54 







5.06 

5.37 



26 

1.42381 



49 




- 



5.13 

6.85 



27 

1.13262 



42 







5.21 

8.47 



28 

1.17520 



34 







5.22 

8.92 



29 

1.09985 



30 







5.20 

8.37 



30 

1.04072 



27 







5.15 

7.39 



31 

.88157 



20 







4.94 

2.92 



32 

.83509 



28 







4.84 

.88 



33 

.60153 



12 







4.77 

-.72 



34 

.53541 



10 







4.44 

-7.64 



35 

.45598 



8 







4.38 

-8.82 



36 

.28576 




1 







3.56 

-25.91 



37 

.12031 



1 





3 

0 

2.02 

-32.82 



38 

.09028 




.8 





3 

0 

2.08 

-30. d3 



39 

.09548 




.6 





2. 

0 

1.50 

-25.23 



40 

.05424 




.4 





2 

0 

1 .66 

-17.20 



41 

.02073 




.2 





1 

5 

1.94 

29.65 



42 

1.49645 

.75749 

20 




5 

6 

5.51 

1.60 

LuDricant 

43 

1.48407 







5 




5.55 

.89 

entrainment 

44 

1.45318 







10 



5.53 

1.25 

direction 

45 

1.40431 







15 



5.50 

1.79 



46 

1.34265 







20 



5.49 

1.96 



47 

1.27083 







25 



5.48 

2.14 



48 

1.19408 







30 



5.52 

1.42 



49 

1.10916 







35 



5.58 

.30 



50 

1.02002 







40 



5.51 

-1.65 



51 

.92775 







45 



5.34 

-4.64 



52 

.10874 

.25249 

16 



0 

5. 

8 

5.89 

1.55 



53 

.14001 

.53866 

20 



0 

15. 

5 

14.5 

-6.45 
















TABLE 7 - COMPARISON BETWEEN 
PRESENT RESULT AND 
EQUATION PROPOSED 


BY HOUPERT 


Case 

H Q (Eq. (39)) 

H Q (Eq. (35)) 

1 

1.52 

2 

1.50 

3 

1.42 

4 

1.38 

10 

1.20 

11 

1.54 

12 

1.46 

13 

1.40 

16 

1.34 

17 

1.34 

18 

1.28 

19 

1.28 

20 

1.28 

21 

1.37: 

24 

1.36 

30 

1.35 

,34 

1.36 

37 

.969 

38 

.968 

39 

.927 





Fig. 1 Lubricant viscosity as a function of 
pressure. Lubricant, sebacate; temper- 
ature, 25 °C. 




(a) Two rigid solids separated by a lubricant film. 



y-Oplane x-Oplane 

(b) Equivalent system of a rigid solid near a plane separated by a lubricant film 

Fig. 2 Contact geometry. 






DIMENSIONLESS MINIMUM FILM 



Fig. 5 Dimensionless minimum film thickness as a function of radius 
ratio. 



(a) a= 4. (b) a = 8. (c) a= 16. 

Fig. 6 Pressure contour for three radius ratios a. 



DIMENSIONLESS FILM PARAMETER, C H min /H min,r 



Fig. 8 Dimensionless film parameter as a function of materials 
parameter. 


DIMENSIONLESS MINIMUM FILM THICKNESS 



Fig. 9 Dimensionless minimum film thickness as a function 
of lubricant entrainment direction. 



Fig. 10 Comparison between minimum film thickness 
formulas. 
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